Sonic analog of gravitational black holes in Bose-Einstein condensates 
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It is shown that, in dilute-gas Bose-Einstein condensates, there exist both dynamically stable and unstable con- 
figurations which, in the hydrodynamic limit, exhibit a behavior resembling that of gravitational black holes. 
The dynamical instabilities involve creation of quasiparticle pairs in positive and negative energy states, as in 
the well-known suggested mechanism for black hole evaporation. We propose a scheme to generate a stable 
sonic black hole in a ring trap. 



03.75.Fi, 04.70.Dy, 04.80.-y 

Many investigations of dilute gas Bose-Einstein conden- 
sates are directed towards experimentally creating nontrivial 
configurations of the semiclassical mean field, or to predict- 
ing the properties of such configurations in the presence of 
quantum fluctuations. Such problems are hardly peculiar to 
condensates, but ultracold dilute gases are so easy to manipu- 
late and control, both experimentally [|l]] and theoretically [^], 
that they may allow us to analyze less amenable systems by 
analogy. As an essay in such an application of condensates, 
in this paper we discuss the theoretical framework and pro- 
pose an experiment to create the analog of a black hole in the 
laboratory and simulate its radiative instabilities. 

The hydrodynamic analog of an event horizon was sug- 
gested originally by Unruh [Qj as a more accessible phe- 
nomenon which might shed some light on the Hawking effect 
[||] (thermal radiation from black holes, stationary insofar as 
backreaction is negligible) and, in particular, on the role of ul- 
trahigh frequencies [0-^]. An event horizon for sound waves 
appears in principle wherever there is a closed surface through 
which a fluid flows inwards at the speed of sound, the flow be- 
ing subsonic on one side of the surface and supersonic on the 
other. There is a close analogy between sound propagation 
on a background hydrodynamic flow, and field propagation 
in a curved spacetime; and although hydrodynamics is only 
a long-wavelength effective theory for physical (super)fluids, 
so also field theory in curved spacetime is to be considered a 
long-wavelength approximation to quantum gravity De- 
termining whether and how sonic black holes radiate sound, 
in a full calculation beyond the hydrodynamic approximation 
or in an actual experiment, can thus offer some suggestions 
about black hole radiance and its sensitivity to high frequency 
physics. 

The basic challenge of our proposal is to keep the trapped 
Bose-Einstein gas sufficiently cold and well isolated to main- 
tain a locally supersonic flow long enough to observe its in- 
trinsic dynamics. Detecting thermal phonons radiating from 
the horizons would obviously be a difficult additional prob- 
lem, since such radiation would be indistinguishable from 
many other possible heating effects. This further difficulty 
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does not arise in our proposal, however, because the black- 
hole radiation we predict is, unlike Hawking radiation, not 
quasistationary, but grows exponentially under appropriate 
conditions. It should therefore be observable in the next gen- 
eration of atom traps. 

A Bose-Einstein condensate is the ground state of a second 
quantized many body Hamiltonian for N interacting bosons 
trapped by an external potential V^, x t (x) [Q] . At zero tempera- 
ture, when the number of atoms is large and the atomic inter- 
actions are sufficiently small, almost all the atoms are in the 
same single-particle quantum state ^(x, t), even if the system 
is slightly perturbed. The evolution of is then given by the 
well-known Gross-Pitaevskii equation 
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where m is the mass of the atoms, a is the scattering 
length, and we normalize to the total number of atoms 
Jd 3 x|^(x,t)| 2 = N. 

Our purposes do not require solving the Gross-Pitaevskii 
equation with some given external potential 14xt(x); our con- 
cern is the propagation of small collective perturbations of the 
condensate, around a background stationary state (x, t) 
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it is only necessary that it be possible, in any external poten- 
tial that can be generated, to create a condensate in this state. 
Many realistic techniques for "quantum state engineering," to 
create designer potentials and bring condensates into specific 
states, have been proposed, and even implemented success- 
fully [[loll; our simulations indicate that currently known tech- 
niques should suffice to generate the condensate states that we 
propose. 

Perturbations about the stationary state 4 r s (x, t) obey the 
Bogoliubov system of two coupled second order differential 
equations. Within the regime of validity of the hydrodynamic 
(Thomas-Fermi) approximation [gj, these two equations for 
the density perturbation q and the phase perturbation qt> in 
terms of the local speed of sound c(x) = — ^/47rap(x), and 



1 



Garay, Anglin, Cirac, and Zoller. Sonic analog of gravitational black holes in Bose-Einstein condensates 



2 



the background stationary velocity v = — Vi3 read 
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Furthermore, low frequency perturbations are essentially just 
waves of (zero) sound. Indeed, the Bogoliubov equations may 
be reduced to a single second order equation for the conden- 
sate phase perturbation <f>. This differential equation has the 
form of a relativistic wave equation d^(yj— gg^ v d v 4>) = 0, 
with g = dot g^, in an effective curved spacetime with the 
metric g^ v being entirely determined by the local speed of 
sound c and the background stationary velocity v. Up to a 
conformal factor, this effective metric has the form 
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This class of metrics can possess event horizons. For in- 
stance, if an effective sink for atoms is generated at the cen- 
ter of a spherical trap (such as by an atom laser out-coupling 
technique JTT|]), and if the radial potential profile is suitably 
arranged, we can produce densities p(r) and flow velocities 
v(x) = — v(r)r/r such that the quantity c 2 — v 2 vanishes at a 
radius r = rh, being negative inside and positive outside. The 
sphere at radius rh is a sonic event horizon completely anal- 
ogous to those appearing in gravitational black holes, in the 
sense that sonic perturbations cannot propagate through this 
surface in the outward direction The physical mecha- 

nism of the sonic black hole is quite simple: inside the hori- 
zon, the background flow speed v is larger than the local speed 
of sound c, and so sound waves are dragged inwards. 

In fact there are two conditions which must hold for this 
dragged sound picture to be accurate. Wavelengths larger 
than the black hole itself will of course not be dragged in, 
but merely diffracted around it. And perturbations must have 
wavelengths A 3> 27r£, 2tt£_/ ^/|1 — v/c\, where £(x) = 
S/[mc(x)] is the local healing length. Otherwise they do not 
behave as sound waves since they lie outside the regime of 
validity of the hydrodynamic approximation. These short- 
wavelength modes must be described by the full Bogoliubov 
equations, which allow signals to propagate faster than the 
local sound speed, and thus permit escape from sonic black 
holes. Even if such an intermediate range of wavelengths does 
exist, the modes outside it may still affect the stability of the 
black hole as discussed below. 

As it stands, this description is incomplete. The conden- 
sate flows continually inwards and therefore at r = there 
must be a sink that takes atoms out of the condensate. Oth- 
erwise, the continuity equation V(/?v) = 0, which must hold 
for stationary configurations, will be violated. We have an- 
alyzed several specific systems which may be suitable the- 
oretical models for future experiments, and have found that 
the qualitative behavior is analogous in all of them. Black 
holes which require atom sinks are both theoretically and ex- 
perimentally more involved, however; moreover, maintaining 
a steady transonic flow into a sink may require either a very 



large condensate or some means of replenishment. We will 
therefore discuss here an alternative configuration which may 
be experimentally more accessible and whose description is 
particularly simple: a condensate in a very thin ring that ef- 
fectively behaves as a periodic one-dimensional system. Un- 
der conditions that we will discuss, the supersonic region in 
a ring may be bounded by two horizons: a black hole hori- 
zon through which phonons cannot exit, and a "white hole" 
horizon through which they cannot enter. 

In a sufficiently tight ring-shaped external potential of ra- 
dius R, motion in radial (r) and axial (z) cylindrical coor- 
dinates is effectively frozen. We can then write the wave 
function as ^(z, r, 9, r) — f(z, r)&(9, r) and normalize $ to 
the number of atoms in the condensate J w d9\<&{9)\ 2 = N, 
where with the azimuthal coordinate 9 we have introduced the 
dimensionless time r = -^jt. The Gross-Pitaevskii equation 
thus becomes effectively one-dimensional: 
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where U = 4vraiVi? 2 / dzdrr\f(z, r)\ 4 and V oxt (0) is the 
dimensionless effective potential (in which we have already 
included the chemical potential) that results from the dimen- 
sional reduction. The stationary solution can then be writ- 
ten as $ s (6*,t) = yj p(9)e l J d6v W an( j the local dimension- 
less angular speed of sound as c{9) = y/Up(ff)/N. Peri- 
odic boundary conditions around the ring require the "wind- 
ing number" w = ^ d9v(9) to be an integer. 

The qualitative behavior of horizons in a ring is well repre- 
sented by the two-parameter family of condensate densities 

N 

p(e) = —(l + bcos9), 

where b E [0, 1]. Continuity, dg(pv) = 0, then determines the 
dimensionless flow-velocity field 
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which depends on w as a third discrete independent parame- 
ter. Requiring that Q s (9, r) be a stationary solution to Gross- 
Pitaevskii equation then determines how the trapping poten- 
tial must be modulated as a function of 9. All the properties 
of the condensate, including whether and where it has sonic 
horizons, and whether or not they are stable, are thus func- 
tions of U, b and w. For instance, if we require that the hori- 
zons be located at 9h = ±7r/2, which imposes the relation 
U = 27rw 2 (l — 6 2 ), then we must have c 2 — v 2 positive for 
9 e (— 7r/2, 7r/2), zero at 9h — ±7r/2, and negative other- 
wise, provided that U < 2irw 2 . The further requirement that 
perturbations on wavelengths shorter than the inner and the 
outer regions are indeed phononic implies U ^> 2n, which in 
turn requires w ^> 1 and 1 3> b ^> 1/w 2 . In fact, detailed 
analysis shows that w > 5 is sufficient. 
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A black hole solution should also be stable over sufficiently 
long time scales in order to be physically realizable. Since sta- 
bility must be checked for perturbations on all wavelengths, 
the full Bogoliubov [||] spectrum must be determined. For 
large black holes within infinite condensates, this Bogoliubov 
problem may be solved using WKB methods that closely re- 
semble those used for solving relativistic field theories in true 
black hole spacetimes [ph. The results are also qualitatively 
similar to those we have found for black holes in finite traps, 
where we have resorted to numerical methods because, in 
these cases, WKB techniques may fail for just those modes 
which threaten to be unstable. 

Our numerical approach for our three-parameter family of 
black/white holes in the ring-shaped condensate has been 
to write the Bogoliubov equations in discrete Fourier space, 
and then truncate the resulting infinite-dimensional eigenvalue 

problem. Writing the wave funtion as $ = $ s + ipe 1 1 d6v ( e \ 
decomposing the perturbation ip in discrete modes 

+ e^e- ine Al^ n {6% 

and substituting into the Gross-Pitaevskii equation, we obtain 
the following equation for the modes w Wj „ and v u n \ 
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In this equation, 
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Eliminating Fourier components above a sufficiently high cut- 
off Q has negligible effect on possible instabilities, which can 
be shown to occur at relatively long wavelengths. The nu- 
merical solution to this eigenvalue equation, together with the 
normalization condition J d9[u^ n Uu;t !n i — u£« n iv,n') = 
Snn'Suiui', provides the allowed frequencies. Real nega- 
tive eigenfrequencies for modes of positive norm are always 
present, which means that black hole configurations are en- 
ergetically unstable, as expected. This feature is inherent in 
supersonic flow, since the speed of sound is also the Landau 
critical velocity. In a sufficiently cold and dilute condensate, 
however, the time scale for dissipation may in principle be 



made very long, and so these energetic instabilities need not 
be problematic Jl2| ] . 

More serious are dynamical instabilities, which occur for 
modes with complex eigenfrequencies and are genuine physi- 
cal phenomena. For sufficently high values of the cutoff (e.g., 
Q > 25 in our calculations), the complex eigenfrequencies 
obtained from the truncated eigenvalue problem become in- 
dependent of the cutoff within the numerical error. The ex- 
istence and rapidity of dynamical instabilities depend sensi- 
tively on (14, b, w). For instance, see Fig. [l] for a contour plot 
of the maximum of the absolute values of the imaginary parts 
of all eigenfrequencies for w — 7, showing that the regions 
of instability are long, thin fingers in the (U, b) plane. Not 
shown in the figure is the important fact that the size of the 
imaginary parts, which gives the rate of the instabilities, in- 
creases starting from zero, quite rapidly with b, although they 
remain small as compared with the real parts. 

The stability diagram of Fig. [j] suggests a strategy for cre- 
ating a sonic black hole from an initial stable state. Within 
the upper subsonic region, the vertical axis 6 = corresponds 
to a homogeneous persistent current in a ring, which can in 
principle be created using different techniques J13]]. Gradu- 
ally changing U and b, it is possible to move from such an 
initial state to a black/white hole state, along a path lying al- 
most entirely within the stable region, and passing only briefly 
through instabilities where they are sufficiently small to cause 
no difficulty. 
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FIG. 1. Stability diagram for winding number w = 7. Solid 
dark-grey areas represent the regions of stability. Smaller plots at 
higher resolution confirm that the unstable "fingers" are actually 
smooth and unbroken. Points on the dashed curve are states with 
horizons at ±7r/2, so that the black/white hole fills half the ring. 

Indeed, we have simulated this process of adiabatic creation 
of a sonic black/white hole by solving numerically (using the 
split operator method) the time-dependent Gross-Pitaevskii 
equation (|l]) that provides the evolution of the condensate 
when the parameters of the trapping potential change so as 
to move the condensate state along various paths in parame- 
ter space. One of these paths is shown in Fig. [I] (light-grey 
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solid line): we start with a current at w = 7, b = 0, and suf- 
ficiently high U; we then increase b adiabatically keeping U 
fixed until an appropriate value is reached; finally, keeping b 
constant, we decrease U adiabatically (which can be physi- 
cally implemented by decreasing the radius of the ring trap), 
until we meet the dashed contour for black holes of comfort- 
able size. Our simulations confirm that the small instabilities 
which briefly appear in the process of creation do not disrupt 
the adiabatic evolution. The final quantum state of the con- 
densate, obtained by this procedure, indeed represents a sta- 
ble black/white hole. We have further checked the stability 
of this final configuration by numerically solving the Gross- 
Pitaevskii equation ([[]) for very long periods of time (as com- 
pared with any characteristic time scale of the condensate) and 
for fixed values of the trap parameters. This evolution reflects 
the fact that no complex frequencies are present, as predicted 
from the mode analysis, and that the final state is stationary. 

Once the black/white hole has been created, one could fur- 
ther change the parameters (Id, b) so as to move between the 
unstable "fingers" into a stable region of higher b (a deeper 
hole); or one could deliberately enter an unstable region. In 
the latter case, the black hole should disappear in an explosion 
of phonons, which may be easy to detect experimentally. Such 
an event might be related to the evaporation process suggested 
for real black holes in the sense that pairs of quasiparticles 
are created near the horizon in both positive and negative en- 
ergy modes. The Hermiticity of the Bogoliubov Hamiltonian 
implies that eigenmodes with complex frequencies appear al- 
ways in dual pairs, whose frequencies are complex conjugate. 
In the language of second quantization, the linearized Hamil- 
tonian for each such pair has the form 

n 

and the only nonvanishing commutators among these opera- 
tors are [A w ,„, A^, ,] = 5 nn i. It is then clear that none of 
these operators is actually a harmonic oscillator creation or 
annihilation operator in the usual sense. However, the linear 
combinations (note that A^, n ^ A]^ n ) 
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and their Hermitian conjugates are true annihilation and cre- 
ation operators, with the standard commutation relations, and 
in terms of these the Bogoliubov Hamiltonian becomes 



state physics [|15|], and optics [16|; but the outstanding recent 
experimental progress in cooling, manipulating and control- 
ling atoms jTo| ] makes Bose-Einstein condensates an espe- 
cially powerful tool for this kind of investigation. We have 
analyzed in detail the case of a condensate in a ring trap, and 
proposed a realistic scheme for adiabatically creating stable 
sonic black/white holes. 

We thank the Austrian Science Foundation and the Eu- 
ropean Union TMR networks ERBFMRX-CT96-0002 and 
ERB-FMRX-CT96-0087. 

Note added. — Further details as well as the study of cigar 
shaped condensates with atom sinks at the center will appear 
inRef. [111. 
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which obviously leads to self-amplifying creation of positive 
and negative frequency pairs. Evaporation through an expo- 
nentially self-amplifying instability is not equivalent, how- 
ever, to the usual kind of Hawking radiation [Bfl; this issue 
will be discussed in detail elsewhere. 

Trapped bosons at ultralow temperature can provide an ana- 
log to a black-hole spacetime. Similar analogs have been pro- 
posed in other contexts, such as superfluid helium Jl4[], solid 



